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We study the localization of protons in the core of neutron stars for ten 
realistic nuclear models that share a common behaviour of nuclear sym- 
metry energy which saturates and eventually decreases at high densities. 
This results in the low proton fraction of beta-stable neutron star matter. 
Protons form a small admixture in the neutron star core, which is localized 
at sufficiently high densities. For every model we calculate the density nioc 
above which the localization effect is present. Our results indicate that 
localization occurs at densities above 0.5 — 1.0 /to~'^. The phase with lo- 
calized protons occupies a spherical shell or a core region inside neutron 
stars which contains significant fraction of all nucleons. Proton localiza- 
tion is of great importance for astrophysical properties of neutron stars 
as it strongly affects transport coefficients of neutron star matter and can 
produce spontaneous magnetization in neutron stars. 



PACS numbers: 21.65. -|-f,97.60.Jd 

Observational parameters of neutron stars shed light on some proper- 
ties of nucleon interactions. These interactions determine the structure and 
properties of neutron stars. The transport and magnetic properties of neu- 
tron stars depend strongly on the structure of the so called liquid core. 
Particularly important is the fraction and structure of the proton compo- 
nent. 

* Presented at the XXIX Mazurian Lakes Conference on Physics, Piaski, Poland, Au- 
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Nuclear symmetry energy plays a crucial role for composition and other 
important properties of dense matter in neutron star. The interaction 
part of the nuclear symmetry energy Eg (n) determines the composition 
of dense nuclear matter Vanishing of the nuclear symmetry energy 
implies proton- neutron separation instability in dense nuclear matter 
Negative values of the symmetry energy result in disappearance of protons 
at high densities. Realistic nuclear matter calculations [^Ej show that the 
symmetry energy becomes negative at high densities, and in consequence 
this leads to the vanishing of protons. 

There are two mechanisms of separation of protons and neutrons in neu- 
tron star matter: a bulk separation of protons and neutrons and localization 
of individual protons in neutron matter. We consider in this paper the lat- 
ter possibility. A bulk separation means that pure neutron matter coexists 
with nucleon matter containing some proton fraction Xc- 

Protons which form admixture tend to be localized in potential wells 
corresponding to neutron matter inhomogeneities created by the protons in 
the neutron medium. To compare the energy of a normal phase of uniform 
density and a phase with localized protons we apply the Wigner-Seitz ap- 
proximation and divide the system into cells, each of them enclosing a single 
localized proton |^ [3] . The neutron background is treated in the Thomas- 
Fermi approximation and the localized proton is described by the Gaussian 
wave function. The neutron density profile is obtained by solving the ap- 
propriate variational equation [^. We calculate the difference of the energy 
per proton in the two phases and look for a critical density for localization 
at which the localized phase becomes the ground state. 

The localization effect is a result of the interaction of protons with small 
density oscillations of the neutron background The protons behave as 
localized polarons which form a periodic lattice at high densities jTj.The aim 
of this paper is to study the proton localization for a number of realistic 
nuclear models with selfconsisted variational method. 

The amount of protons present in the neutron star matter, which is 
charge neutral and /^-stable, is crucial for the cooling rate of neutron stars 
and also plays an important role for magnetic and transport properties of 
neutron star matter. Nuclear models do not uniquely predict the proton 
fraction of the neutron star matter at high densities. This controversy is 
discussed in details in Ref. ^ |H| where the discrepancy of the proton frac- 
tion in various models is shown to reflect the uncertainty of the nuclear 
symmetry energy at high densities. In this paper we consider a class of 
nuclear interaction models for which the proton fraction is of the order of a 
few (for five of them of above ten) percent and decreases at high densities. 
For the calculations we have chosen ten realistic nuclear interaction models. 
These are interactions derived by Myers and Swiatecki |^ (MS), the Skyrme 



A.Szmaglinski printed on February 5, 2008 



3 



model with parameters (SI', SII', SIII', SL, Ska, SKM, SGII, RATP, T6) 
from Ref . jlUl lllj , the Friedman and Pandharipande interactions ^| (as 
parametrized by Ravenhall in Ref-PUj) (FPR), three models, UV14+TNI, 
AV14+UVII and UV14+UVII, from Ref.f^ by Wiringa et al. and four mod- 
ern APR models , which provides a fit to the nucleon-nucleon scattering 
data in the Nijmegen data base. Fig.l shows the proton fraction of beta- 
stable and charge neutral matter containing neutrons, protons, electrons 
and muons for all the interactions. 
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Fig. 1. Proton fraction of neutron star matter for indicated interaction models. 

In calculactions we apply selfconsisted method described in our paper . 
The threshold densities (in fm~^) as well as rms proton radii (in fm), when 
the localization effect becomes energetically favorable are listed in Table 1. 
We also present the minimum masses of neutron stars Mmirn when central 
density is equal to the threshold density of localization of protons. For every 
model we calculated also the maximum mass of neutron stars Mmax (in solar 
units). The localization effect is present above the threshold barion number 
density nioc- For some models protons vanish from the system above some 
density n^, so the localization effect is present in the spherical shape region 
of neutron stars in the density range {nioc',ny) Fig. 2. For Ska parametriza- 
tion of the Skyrme model and the A18-|-5v model the threshold density for 
proton localization is also the central density of predicted maximum mass 
configuration of neutron star, so there is no localization of proton inside the 
observed neutron stars in these two cases. 

Localization of protons influences the global properties of neutron star, 
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SO it is of importance to calculate the fraction of nucleons inside the region 
with localized protons in the neutron star. We show this fraction in Fig. 3. 

Results of our calculations for nuclear interactions we use indicate that 
the proton impurity in neutron star matter becomes localized at densities 
above 0.5 — 1.0/m~^. This has important consequences for neutron stars 
as densities in this range correspond to the inner core of neutron stars with 
masses exceeding one solar mass, M > IMq. 

TABLE 1 

For considered models we present the threshold densities for localization and correspond- 
ing rms radii of proton wave function calculated using selfconsistent variational method, 
the minimum mass of neutron star Mmin (in solar units) which has central density equal 
to the threshold density of localization and the predicted maximum mass od neutron star. 



Potential 


nioc 


Ttloc 
Tip 




-^max 


MS 


1.030 


0.905 


2.019 


2.038 


sr 


0.351 


1.519 


0.876 


2.242 


sir 


0.361 


1.688 


0.829 


1.966 


siir 


0.337 


1.552 


0.898 


2.266 


SL 


0.964 


1.384 


1.480 


1.631 


Ska 


1.016 


0.804 


2.241 


2.241 


SKM 


0.979 


1.330 


1.570 


1.677 


SGII 


0.899 


1.367 


1.432 


1.661 


RATP 


0.709 


1.397 


1.225 


1.726 


T6 


0.530 


1.985 


0.600 


1.429 


FPR 


0.721 


1.262 


1.435 


1.800 


UV14+TNI 


0.731 


1.209 


1.489 


1.827 


AV14+UVII 


0.789 


0.971 


1.793 


2.123 


UV14+UVII 


0.766 


0.913 


1.986 


2.207 


A18 


1.493 


1.136 


1.647 


1.673 


A18+5v 


1.627 


0.915 


1.805 


1.805 


A18+UIX 


0.645 


0.911 


2.095 


2.386 


AlS+^v+UIX* 


0.819 


0.878 


2.052 


2.213 



According to the Ref. significant number of measured neutron star 
masses are inside the range (Mmin; -^max)- It indicates that localization of 
protons is an universal state of dense nuclear matter in neutron stars. 

Strongly asymmetric nuclear matter is unstable with respect to proton 
localization jl5l U|. A uniform proton distribution and a periodic proton 
arrangement result in very different properties The presence of the 
localized protons inside neutron star cores would have profound astrophys- 
ical consequences. In particular, the cooling proceeds in a quite different 
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Fig. 2. The spherical shell of neutron star, where the localization effect occurs. The 
region corresponds to density above the threshold and below the density for which 
the protons vanish. Neutron star density profiles are labelled by logarithm of the 
central density. 
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Fig. 3. The fraction of nucleons in the region with localized protons of the neutron 
star versus neutron star mass. 



way. Recent analysis jl6| shows that the presence of such localized proton 
phase results in more satisfactory fits of temperatures of observed neutron 
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stars. Also, spin ordering of localized protons could strongly affect magnetic 
properties of the system ^1 E| • The spin ordered phase can contribute sig- 
nificantly to the observed magnetic moments of neutron stars jl8l I19j . 

This research is partially supported by the Polish State Committee for 
Scientific Research, grants 2P03B 110 24 and PBZ-KBN-054/P03/2001. 
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